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ABSTRACT 


The  Refractivity  Structure  Matching  Algorithm  (RSMA)  is  an  automated  aid  for  matching 
refractivity  stmctures  in  adjacent  profiles  which  result  from  common  layers.  The  RSMA  is 
specifically  designed  to  prepare  multiple  refractivity  profiles  for  entry  into  the  Naval  Ocean 
Systems  Center's  Radio  Physical  Optics  (RPO)  program.  A  full  description  of  the  RSMA,  a  test 
case  demonstrating  its  capabilities,  and  a  computer  listing,  are  provided  in  this  document. 
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SPECIFICATION  FOR  A  REFRACTTVITY  STRUCTURE  MATCHING  ALGORITHM 


1.  INTRODUCTION 

Refractivity  layers  in  the  lower  atmosphere  result  from  rapid  vertical  changes  in  air 
temperature  and/or  humidity.  Variations  in  refractive  structures  over  altitude  and  range  may  sig¬ 
nificantly  affect  electromagnetic  propagation.  In  order  to  represent  the  evolution  of  refractivity 
layers  with  range,  the  Refractivity  Structure  Matching  Algorithm  (RSMA)  has  been  developed. 
TTie  RSMA  is  specifically  designed  to  prepare  multiple  profiles  of  M-units  versus  height,  as  a 
function  of  range,  for  entry  into  the  NOSC  range-dependent  Radio  Physical  Optics  (RPO)  model 
(NOSC,  1991). 

The  methodology  used  in  the  RSMA  to  match  up  refractive  structures  in  adjacent  profiles 
that  result  from  common  layers  follows  procedures  developed  at  the  John  Hopkins  AFL 
(Konstanzer  and  Dockery,  1989;  Konstanzer,  1989).  Prior  to  matching,  each  point  on  each 
refractivity  profile  is  categorized  as  being  one  of  31  structural  types,  based  on  the  profile 
M-gradients  immediately  below  and  above  the  point,  and  the  magnitude  and  direction  of  the 
change  in  these  gradients  across  the  point.  Structure  ty^jes  in  adjacent  profiles  are  ne.\t  matched 
up  in  a  specific  order  based  on  a  hierarchy,  such  that  the  structures  considered  most  likely  to 
affect  propagation  are  matched  first.  Data  points  not  initially  matched  by  the  hierarchal  pairings 
are  matched  by  the  addition  of  points  to  the  profiles  by  interpolation.  This  secondary  piPtching 
ensures  that  two  important  RPO  data  requirements  are  met,  namely,  1)  every  profile  l.as  the  same 
number  of  points  and  2),  the  nth  point  in  any  one  profile  is  matched  to  the  nth  points  in  the  other 
profiles. 


2. 


ALGORITHM  FUNCTIONS 


2.1  Inputs,  Outputs  and  Limits 

The  following  parameters  and  refractivity  field  data  are  required  by  the  RSMA  and  need 
to  be  specified  by  the  operator: 

1)  Number  of  M-unit  versus  height  profiles  -  2  to  11 

2)  Number  of  data  points  for  each  profile  -  3  to  30 

3)  Couplets  of  M-units  and  height  for  each  profile,  corresponding  to  profile  data  points. 

4)  Height  unit  -  ’  m'  or  ’ft' 

5)  Range  for  each  profile  -  >=  0. 

6)  Range  unit  -  ’  km’,  ’mi’  or  ’nmi’ 

The  first  profile  specified  must  be  for  a  range  of  zero;  the  range  of  subsequent  profiles 
must  be  greater  than  the  range  of  the  previously  entereJ  profile.  The  RSMA  places  no  numerical 
limits  on  either  height  or  M-unit  values;  however,  height  array  values  must  be  specified  in 
ascending  order  and  all  height  profiles  must  have  the  same  value  for  the  first  (and  last)  array 
elements. 

After  RSM.A  execution,  profiles  defined  by  height  and  M-unit  data  couplets  are  specified, 
along  with  the  number  of  profile  data  points  (i.e.,  profile  levels).  As  required  by  the  RPO  pro¬ 
gram,  all  RSMA  output  profiles  have  the  same  number  of  points  (levels),  with  the  nth  point  in 
any  one  profile  matched  to  the  nth  point  in  the  other  profiles.  In  general,  the  number  of  profile 
levels  determined  by  the  RSMA  depends  upon  the  number  of  profiles,  the  number  of  profile 
pioints,  and  the  number  of  common  refractivity  features  matched  between  profiles.  A  priori,  this 
last  element  is  not  known.  Operationally,  to  have  a  decent  likelihood  in  meeting  the  current  RPO 


requirement  of  no  more  than  51  levels,  the  maximum  possible  number  of  levels  after  RSNIA 
execution,  nlvmx.  given  by 

nlvmx  =  S  npt(i)  -  (nprf*2)  +  2  , 

where  npt(i)  is  the  number  of  points  in  the  ith  profile  and  nprf  is  the  number  of  profiles  , 
should  not  exceed  the  RPO  limit  by  more  than  50%  (nlvmx  =  ~  75)  prior  to  RSMA  execution. 

Output  profiles  correspond  to  those  originally  specified  by  the  operator  in  regards  to 
number,  range  and  units.  A  sample  RSMA  output  is  presented  in  the  Appendix;  these  same  data 
are  depicted  graphically  in  Section  3.  Graphics  and  profile  interpolation  to  intermediate  ranges 
are  not  included  with  the  RSMA;  these  functions  are  provided  by  the  RPO  program. 

2.2  Data  Point  Categorization 

Refractivity  layers  in  two  adjacent  profiles  may  be  matched  provided  the  layers  have  the 
same,  or  similar,  structure  categories.  In  the  RSMA,  categoiies  are  assigned  to  each  data  point 
of  each  refractivity  profile,  based  on  profile  M-gradients  below  and  above  the  point,  and  the 
change  of  these  gradients  across  the  point.  For  the  nth  point  of  the  ith  profile,  the  M-gradient 
(in  M  units/km  )  is  defined  as 
dM/dz(i.l)  =  0. 
and 

dM/dz(in)  =  1000.  *  (  M(i,n)-M(i,n-1)  /  z(i.n)-z(i.n-l)  ) 
where  M  is  the  modified  refractivity  in  M-units,  and  z  is  the  height  in  meters. 

The  M-gradient  defines  the  four  basic  tv-pes  of  refraction: 

Trapping  -  TRP  <  0  M-units/km 

Sujjerrefractive  -  SUP  >=  0  to  <  78  M-units/km 
Standard  -  STD  >=  78  to  <  157  M-units/km 

Subrefractive  -  SUB  >=  157  M-units/km 


The  change  in  the  M-gradieni  across  the  nth  point  of  the  ith  profile  (in  M-units/Tcm)  is  defined 
as 

A  (dM/dz)(iji)  =  0.  for  n  =  1  and  nmx 

where  nmx  is  the  last  (highest  indexed)  data  point,  and 

A  (,dM/dz)(i.n)  =  dM/dz(i.n+l)  -  (iM/dz(i,n)  for  n  =  [2,nmx-l] 

Table  1  presents  the  31  data  point  structure  upes.  These  categories  are  given  by  integers 
over  the  range  [-15.15].  Negative  (Positive)  categories  indicate  a  decrease  (increase)  in  M- 
gradient  acro.«s  a  data  point.  Positive  and  negative  categories  of  like  absolute  value  have  been 
specified  as  "mirror"  images  of  one  another.  The  categorization  of  profile  data  points  is  based 
on  a  hierarchy,  chosen  to  characterize  which  stmctures  are  most  likely,  and  least  likely,  to  affect 
propagation.  The  structure  type  T5’  (‘-15’),  which  defines  a  M-gradient  change  from  trapping 
to  subrefractive  (subrefractive  to  trapping)  across  a  data  point,  is  considered  to  be  of  most 
importance.  On  the  other  hand,  the  structure  type  ’O’,  which  signifies  an  absolute  change  in  M- 
gradient  of  less  than  39  M-units/km  across  a  data  point,  is  considered  to  be  the  least  important. 
The  first  and  last  points  of  all  profiles  are  assigned  the  structure  ty'pe  ’15’;  this  ensure.:  that  the 
bottoms  and  tops  of  adjacent  profiles  are  always  matched  up. 

2.3  Matching 
2.3.1  Primary 

Refractivity  structures  in  adjacent  profiles  may  be  matched  to  each  other  provided  that 
they  are  of  like  or  similar  categories,  and  fail  within  a  certain  altitude  window.  A  hierarchal 
classification  of  paired  structure  types,  given  in  Table  2.  is  used  to  ensure  that  those  refractivity 
structures  deemed  most  important  are  first  matched.  This  hierarchy  defines  the  order  in  which 
structures  in  one  profile  (A)  are  to  be  matched  to  structures  in  the  adjacent  profile  (B).  The  first 
30  pairs  jf  (profile  A.  profile  B )  categories  match  exact  structure  ty-pes  between  the  two  profiles. 
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Table  1.  Profile  structure  t^pe  categories  for  the  nth  point  of  the  ith  profile,  based  on  M- 
gradients  below  and  above  the  data  point,  and  the  change  in  the  M-gradient  across  the 
point.  M-gradients  are  given  as  TRP,  SUP,  STD  and  SUB. 


dM/dzUdi) 

dM/dz(iji+I) 

A  (<lM/dz)(i,n) 

M-units/km 

(below) 

(above) 

< -78 

(-78,-39]  (-39,39) 

[39,78) 

>78 

SUB 

TRP 

-15 

STD 

TRP 

-14 

SUB 

SUP 

-13 

SUP 

TRP 

-12 

STD 

SLB 

-11 

SLB 

STD 

-10 

TRP 

TRP 

-  9 

SUB 

SLB 

-  8 

SUP 

TRP 

-  7 

STD 

SLB 

-  6 

SLB 

STD 

-  5 

TRP 

TRP 

-  4 

SUP 

SUP 

-  3 

STD 

STD 

SLB 

SUB 

-  1 

0 

SLB 

SUB 

1 

STD 

STD 

n 

SUP 

SUP 

3 

TRP 

TRP 

4 

STD 

SUB 

5 

SLB 

STD 

6 

TRP 

SUP 

7 

SLB 

SUB 

8 

TRP 

TRP 

9 

STD 

SUB 

10 

SUP 

STD 

11 

TRP 

SUP 

12 

SLB 

SLB 

13 

TRP 

STD 

14 

TRP 

SLB 

15 
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Table  2.  Hierarchy  of  the  1 14  pairs  of  structure  t\-pes  for  matching.  Note  that  the  first  15  match 
combinations  have  2  permutations,  and  the  last  21  have  4  permutations. 


MATCH  COMBINATION 

1 

PERMUTATIONS 

2  3 

4 

1 

-15  -15 

15  15 

2 

-14  -14 

14  14 

3 

-13  -13 

13  13 

4 

-12  -12 

12  12 

5 

-11  -11 

11  11 

6 

-10  -10 

10  10 

7 

-9  -9 

9  9 

8 

-8  -8 

8  8 

9 

-7  -7 

7  7 

10 

-6  -6 

6  6 

1 1 

-5  -5 

5  5 

12 

-4  -4 

4  4 

13 

-3  -3 

3  3 

14 

-2  -2 

2  2 

15 

-1  -1 

1  1 

16 

-15  -14 

15  14 

-14  -15 

14 

15 

17 

-15  -13 

15  13 

-13  -15 

13 

15 

18 

-14  -13 

14  13 

-13  -14 

13 

14 

19 

-14  -12 

14  12 

-12  -14 

12 

14 

20 

-14  -11 

14  11 

-11  -14 

11 

14 

21 

-13  -10 

13  10 

-10  -13 

10 

13 

-13  -11 

13  11 

-11  -13 

11 

13 

23 

-15  -12 

15  12 

-12  -15 

12 

15 

24 

-15  -10 

15  10 

-10  -15 

10 

15 

25 

-15  -11 

15  11 

-11  -15 

1 1 

15 

26 

-14  -10 

14  10 

-10  -14 

10 

14 

27 

-13  -12 

13  12 

-12  -13 

12 

13 

28 

-12  -11 

12  11 

-11  -12 

11 

12 

29 

-11  -10 

11  10 

-10  -11 

10 

11 

30 

-12  -9 

12  9 

-9  -12 

9 

12 

31 

-10  -8 

10  8 

-8  -10 

8 

10 

32 

-12  -7 

12  7 

-7  -12 

7 

12 

33 

-11  -6 

11  6 

-6  -11 

6 

11 

34 

-10  -5 

10  5 

-5  -10 

5 

10 

35 

-9  -4 

9  4 

-4  -9 

4 

9 

36 

-8  -1 

8  1 

-1  -8 

1 

8 
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The  remaining  84  pairs,  which  match  structures  of  similar  categories,  have  been  ranked  according 
to  a  subjective  judgment  as  to  how  individual  pairings  likely  affect  propagation.  In  a  general 
way,  the  114  hierarchal  pairings  of  Table  2  may  be  construed  as  being  36  principal  match 
combinations,  of  which  the  first  15  contain  two  permutations,  and  the  remaining  21,  four 
permutations  each.  The  numerous  pairs  of  structure  match  categories  which  do  not  appear  in 
Table  2,  such  as  (-7.-4)  or  (7,6),  are  considered  minor  pairings,  and  are  not  used  by  the  RSMA. 

Given  two  adjacent  profiles  (A  and  B),  the  RSMA’s  primary  matching  scheme  begins  by 
searching,  from  the  bottom  and  working  upward,  profile  A  data  points  for  the  stmcture  tvpe 
which  corresponds  exactly  to  the  appropriate  (i.e.,  profile  A)  value  of  the  first  hierarchal  match 
pairing.  If  the  appropriate  structure  type  carmot  be  found  in  profile  A,  the  next  pair  of  match 
categories  are  considered,  and  a  similar  search  is  made  for  the  profile  A  structure  type 
corresponding  to  that  hierarchal  pairing.  This  search  process  terminates  when  all  114  pairs  of 
match  categories  have  been  tried  and  is  only  interrupted  when  the  search  successfully  finds  the 
required  profile  A  structure  needed  for  matching. 

Given  an  appropriate  profile  A  match  category,  an  al'itude  w'indow  is  next  determined  for 
matching.  This  window  requires  that,  for  matching,  structures  in  adjacent  profiles  be  within 
reasonable  physical  (i.e..  vertical)  bounds.  The  size  of  this  window  is  variable,  and  is  a  function 
of  both  the  range  between  profiles  and  the  altitude  of  the  structure  being  matched.  Specifically, 
the  range  dependence  of  the  window  (RDEP)  is  given  by 
RDEP  =  RFCT  *  A  range 

where  RFCT  is  a  constant  corresponding  to  an  altitude  over  range  slope  of  1/1000,  and 
A  range  is  the  difference  in  range  between  two  adjacent  profiles. 

For  A  range  in  km  (nmi)  and  RDEP  in  m  (ft),  RFCT  =  1.00  (6.08). 
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The  altitude  dependence  of  the  window  (ZDEP)  is  given  by 
ZDEP(ft)  =  200.  for  z(i.n)  <  100  ft,  and 

ZDEP(ft)  =  100.  *  logj^z(iji)  for  z(i,n)  >=  100  ft  . 
where  z(i,n)  is  the  altitude  of  the  profile  A  structure  to  be  matched.  ZDEP  values  in  meters 
correspond  exactly  to  these  fomiulae. 

The  altitude  window  for  structure  t^'pe  matching  between  adjacent  profiles  is  given  suuply 
as 

^V'NDW  =  +  (RDEP  +  ZDEP) 

In  general,  the  greater  the  distance  between  profiles  and.  the  higher  the  stmcture  to  be  matched, 
the  larger  the  altitude  window.  As  an  example  of  a  WNDW  calculation,  consider  z(i.n)  =  1000 
ft  and  A  range  =  25  nmi.  Here.  W'NDW  would  be  +  452  ft.  or  from  an  altitude  of  548  to 
1452  ft. 

Once  a  window  has  been  established,  the  RSMA  attempts  to  find  a  profile  B  structure  type 
within  that  window  which  corresponds  to  the  appropriate  category  match  pair.  If  the  search  is 
successful,  the  profile  A  and  profile  B  structure  u^es  (data  points)  are  considered  matched.  If 
more  than  one  profile  B  structure  type  within  the  window  is  a  match,  then  the  profile  B  strucmre 
tv^pe  (data  point)  nearest  in  elevation  to  the  profile  A  data  point  (z(i,n))  is  chosen.  If  a  profile 
B  category  match  cannot  be  found  within  the  window,  the  next  profile  A  data  point  (at  height 
z(i.n+l))  is  considered  for  matching.  Already  matched  profile  stmctures  place  an  important 
limitation  on  subsequent  matching.  Specifically,  given  an  existing  match  between  the  nth  data 
point  of  the  ith  profile  and  the  mth  data  point  of  the  i+l  profile,  the  RSMA  does  not  permit  an 
ith  profile  data  point  above  (below)  the  index  n  to  be  matched  with  a  i+1  profile  data  point  below 
(above)  the  index  m. 
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Once  the  primary  matching  procedure  has  been  completed  for  the  first  two  profiles,  it  is 
repeated  again  for  the  second  and  third  profiles,  and  so  forth,  until  matching  has  been  attempted 
between  all  adjacent  profiles.  Upon  conclusion  of  this  part  of  the  RSMA,  it  is  possible  to  have 
defined  refractivity  structures  common  to  several  profiles  and  extending  over  considerable  hori¬ 
zontal  extent,  or  structures  only  common  to  two  profiles  and  of  limited  horizontal  range. 

2.3.2  Secondary 

The  RSMA  secondary  matching  procedure  matches  all  profile  data  points  which  remain 
unmatched  initially.  Such  data  points  either  correspond  to  minor  refractivity  structures  or.  signifi¬ 
cant  structures  which  fade  rapidly  over  range  (i.e.,  from  one  profile  to  the  next).  Secondary^ 
matching  is  essentially  an  inteipolative  process,  in  which  an  unmatched  data  point  in  a  particular 
profile  is  matched  to  new'  data  points  in  all  the  other  profiles.  On  conclusion  of  this  process, 
every  profile  will  have  the  same  number  of  data  points,  with  the  nth  data  point  in  any  one  profile 
matched  to  the  nth  points  in  the  other  profiles. 

Secondary  matching  is  done  in  two  directional  sweeps,  first  forw  ard,  then  backward,  starting 
with  the  first  two  profiles  (  A  and  B,  A  at  zero  range).  Assume  that  there  exists  an  unmatched 
profile  A  data  point  (structure)  at  a  height  ZI  located  between  previously  matched  (b>'  the 
primary  matching  procedure)  profile  A  data  points  at  heights  ZBl  (below)  and  ZTI  (above). 
Further  assume  that  the  profile  A  points  ZBl  and  ZTI  have  been  previously  matched  (by  the 
primary  matching  routine)  to  profile  B  data  points  ZB2  and  ZT2,  and  that  there  exists  no  other 
primary  matched  stmctures  between  ZBI-ZB2  and  ZT1-ZT2.  Then,  a  profile  B  height  Z2, 
corresponding  to  a  new  match  with  the  profile  A  data  point  at  Zl,  is  obtained  by 
Z2  =  [  (  (Z1-ZB1)/(ZT1-ZB1)  1  *  (ZT2-ZB2)  ]  +  ZBl  . 

This  formula  is  applied  to  all  profile  A  unmatched  data  points.  On  conclusion,  the  new 
interpolated  profile  B  altitude  values  are  inserted  (by  means  of  soning)  into  the  profile  B  height 
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array  at  the  appropriate  locations.  Provided  that  there  are  at  least  three  profiles,  the  secondary 
matching  procedure  is  next  repeated,  with  the  newly  augmented  second  profile  as  profile  A  and 
the  third  profile  as  profile  B,  and  so  forth,  until  the  final  profile  is  reached.  As  tlie  procedure 
steps  forward,  the  number  of  unmatched  data  points  at  any  particular  ith  profile  A  will  be  greater 
than  the  number  for  the  previous  (ith-1)  profile  A,  provided  that  the  ith  profile  itself  has  data 
points  not  initially  matched  by  the  primary  matching  routine.  Once  the  last  profile  (at  the 
maximum  range)  has  been  reached,  the  secondary  matching  process  is  reversed,  the  last  profile 
becomes  profile  A  and  the  next  to  la.st,  profile  B.  Upon  completion  of  matching  between  these 
two  profiles,  backward  matching  is  repeated,  with  the  newly  augmented  profile  B  as  the  new 
profile  A  and  the  third  to  last  profile  as  profile  B,  and  so  forth,  until  the  first  profile  (  at  zero 
range)  has  again  been  reached.  At  this  point,  aU  profiles  will  have  the  same  number  of  data 
points  (i.e.,  levels),  and  the  nth  point  of  any  profile  will  correspond  to  the  nth  points  of  the  other 
profiles.  Redundant  data  points  (levels)  are  next  removed  from  the  profile  arrays.  Finally, 
refractivity  M-values  are  calculated  for  all  newly  interpolated  height  array  values,  by  means  of 
vertical  interpolation  of  existing  M-values  to  these  heights. 

3.  TEST  CASE 

To  test  the  RSMA,  high  resolution  refractivity  field  data  (M-profiles),  corresponding  to 
a  Pt.  Loma  to  Guadalupe  Island  transect  of  12  March  1948,  were  utilized.  This  data,  shown  in 
Figure  1,  consist  of  five  M-units  versus  height  profiles,  at  ranges  of  0,  80,  120,  160  and  200  nmi. 
All  profiles  extend  from  sea  level  to  a  height  of  3500  ft.  The  number  of  profile  data  points  is 
28  for  each  of  the  fnst  three  profiles,  and  24  each  for  the  tw'o  profiles  at  the  greatest  ranges. 
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Figure  1.  Refractivity  field  data,  for  a  Pt.  Loma  to  Guadalupe  Is.  transect  of  12  March  1948. 
M-profiles  are  at  ranges  of  0,  80,  120.  160  and  200  nmi.  Profile  values  are  given  by  the 
projection  of  the  M-scale  onto  the  horizontal  at  the  desired  range. 


11 


Figure  2  depicts  results  from  the  RSMA  after  primary-  matching.  Here,  line  segments  denote 
matched  structures;  note  that  there  are  many  profile  data  points  which  were  not  matched.  TTie 
most  prominent  refractivity  feature  matched  is  a  trapping  layer  sloping  upward  from  about  600 
ft  at  zero  range  to  around  1600  ft  at  200  nmi  range.  Other  refractivity  structures  are  observed 
not  to  extend  over  all  ranges.  The  strong  trapping  layer  near  2700  ft  at  range  200  nmi  is  seen 
to  fade  rapidly  over  range,  and  can  only  be  matched  40  nmi  inward.  As  required  by  the  RSM.'^, 
the  bottoms  and  tops  of  all  profiles  are  matched. 

V^Tiile  the  overall  matching  appears  quite  favorable,  the  matching  of  the  1 6th  data  point  of 
the  second  profile  (at  a  height  of  1560  ft)  with  the  15th  point  of  the  third  profile  (at  a  height  of 
1500  ft),  both  of  structure  category  '12’,  likely  would  not  have  been  done  by  an  expert  analyst. 
More  likely,  such  an  expert  would  have  matched  the  10th  point  of  the  second  profile  (at  a  height 
of  1050  ft,  and  structure  category  ’7’)  with  the  15th  point  of  the  thu-d  profile.  While  a  (7,12) 
category  pairing  does  appear  in  the  RSMA  primary  matching  hierarchy,  it  is  ranked  well  below 
the  (12,12)  hierarchal  pairing. 

Figure  3  depicts  results  from  the  RSMA  after  both  primary  and  secondary-  matching.  Here, 
all  profiles  have  the  same  number  of  data  points  (levels),  98,  and  the  nth  point  in  any  one  profile 
is  matched  to  the  nth  points  of  the  other  profiles.  Figure  3  shows  the  spacing  between  matched 
profile  levels  is  quite  variable,  owing  to  the  spatial  distribution  of  unmatched  data  points  after 
primary  matching.  Tabular  height  and  refractivity  values  corresponding  to  the  matched  profile 
levels  of  Figure  3  can  be  found  in  the  Appendix;  they  are  output  from  the  RSMA.  These 
particular  height  and  refractivity  arrays  (of  98  elements)  would  not  be  suitable  for  direct  entry 
into  the  RPO  program,  since  that  program  limits  the  number  of  levels  to  no  more  than  51.  To 
meet  such  a  limit,  it  would  have  been  necessary  to  "thin"  (i.e.,  decrease  the  number  of  data  points 
of)  the  original  M -profiles.  le.s.sen  the  number  of  profiles,  or  do  a  combination  of  both. 
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Figure  3.  Matched  profile  structures,  indicated  by  line  segments,  after  RSMA  primary  and 
secondary  matching,  for  the  Pt.  Loma  to  Guadalupe  Is.  transect  of  12  March  1948. 


14 


Interpolation  of  processed  M-profiles  to  intermediate  ranges  is  not  a  function  of  the  RSMA, 
since  this  tasK  is  performed  by  the  RPO  program.  Nonetheless,  during  the  developmental  stage 
of  the  RSMA,  linear  interpolation  of  M-profdes  in  range  between  points  in  adjacent  matched 
profiles  was  performed.  Figure  4  depicts  profiles  at  intermediate  ranges,  at  increments  of  8  nmi. 
Such  a  graphical  display  allows  the  viewer  to  clearly  see  the  evolution  of  refractive  layers  over 
range,  and  points  out  the  heterogeneous,  range-dependent  nature  of  large-scale  refractivity  fields. 


320  h-UNITS  420 


RANGE  (NMI ) 

Figure  4.  Origmal  and  intermediate  (derived)  M-profiles,  for  the  Pt.  Loma  to  Guadalupe  Is. 
transect  of  12  .March  1^48.  Intermediate  profiles  are  interpolated,  using  the  final  matched  profile 
structures  from  the  RSMA  (  Fig.  3). 
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APPENDIX  -  RSMA  COMPUTER  LISTLNG 


The  RSMA  is  written  in  ANSI  Fortran  77. 

The  output  from  the  DEMO  Program  follows  the  computer  code. 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


RSMA 

REFRACTIVITY  STRUCTURE  MATCHING  ALGORITHM 


cccccccccccccccccccc 

MAIN  PROGRAM  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

INPUTS 

VARIABLE 

DESCRIPTION 

(VALID 

RANGE, UNITS) 

c 

NPRF 

NO.  OF  PROFILES 

2-11 

c 

NPT(I) 

NO.  OF  POINTS,  ITH  PROFILE 

3-30 

c 

Z(I,NPT(I)  ) 

HEIGHT  ARRAY,  ITH  PROFILE 

c 

ZUNIT 

HEIGHT  UNIT 

'  M'  OR  'FT' 

c 

R(I,NPT(I) ) 

REFRACTIVITY  ARRAY,  ITH  PROFILE 

M-UNITS 

c 

RANGE ( I ) 

RANGE,  ITH  PROFILE 

>=  0. 

c 

RUN  IT 

RANGE  UNIT 

'  KM' , ' 

MI'  OR  'NMI' 

u 

c 

OUTPUTS 

VARIABLE 

DESCRIPTION 

(VALID 

RANGE, UNITS) 

c 

IPMX 

NO.  OF  PROFILE  LEVELS 

c 

AFTER  RIA  PROGRAM  EXECUTION 

c 

ZI(I,I?MX) 

HEIGHT  ARRAY,  ITH  PROFILE 

c 

RI(I,I?MX) 

INTERPOLATED  REFRACTIVITY 

M-UNITS 

c 

ARRAY,  ITH  PROFILE 

c 

C  THE  RSMJ\  IS  DESIGNED  TO  PREPARE  MULTIPLE  M-UNIT  PROFILES  (GIVEN  AS 
C  A  FUNCTION  OF  RANGE)  FOR  ENTRY  INTO  THE  NOSC  RADIO  PHYSICAL  OPTICS 
C  (RPO)  PROGRAM.  THE  RPO  REQUIRES  ALL  PROFILES  TO  HAVE  THE  NUMBER 

C  OF  LEVELS,  AND  LIMITS  THE  NUMBER  OF  LEVELS  TO  NO  MORE  THAN  51.  IN 
C  GENERAL,  THE  NUMBER  OF  PROFILE  LENTILS  DETERMINED  BY  THE  RSMA  DEPENDS 
C  UPON  THE  NUMBER  OF  PROFILES,  THE  NUMBER  OF  PROFILE  POINTS,  AND  THE 
C  NUMBER  OF  COMMON  REFRACTIVITY  FEATURES  MATCHED  BETWEEN  PROFILES.  A 
C  PRIORI,  THIS  LAST  ELEMENT  IS  NOT  KNOWN.  OPERATIONALLY,  FOR  A  DECENT 
C  LIKELIHOOD  IN  MEETING  THE  RPO  REQUIREMENT  OF  NO  MORE  THAN  51  LEVELS, 
C  THE  MAXIMUM  POSSIBLE  NUMBER  OF  LEVELS  AFTER  RSMA  EXECUTION,  NLVMX, 

C  GIVEN  BY  I  NPT(I)  -  (NPRF*2)  +  2  ,\SHOULD  NOT  EXCEED  THIS  LIMIT 
C  BY  MORE  THAN  50%  (NLVMX  =~75)  PRIOR  TO  RSMA  EXECUTION. 

C  THE  FOLLOWING  PROGRAM  IS  A  DEMONSTRATION  DRIVER  FOR  THE  RSMA  SUB- 
C  ROUTINE.  THE  INPUT  DATA  SET  CONSISTS  OF  HIGH  RESOLUTION  REFRACTIVITY 
C  PROFILES  FROM  GUADALUPE  ISLAND.  SI.NCE  THE  PROFILES  ARE  GIVEN  IN 
C  N-UNITS,  THEY  WERE  CONVERTED  TO  M-UNITS  PRIOR  TO  CALL  TO  THE  RSMA 
C  SUBROUTINE.  FOR  THIS  PARTICULAR  DEMONSTRATION  DATA  SET,  THE  ACTUAL 
C  NUMBER  OF  PROFILE  LEVELS  DETERMINED  BY  THE  RSMA  (IPMX  =  98)  SIGNIFI- 
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CANTLY  EXCEEDS  THE  RPO  LIMIT  OF  51  LEVELS. 

START  DEMO  PROGRAM 

PROGRAM  RFSTMT 

CHARACTER* 2  ZUNIT 
CHARACTER* 3  RUNIT 

ENTER  PROFILE  PARAMETERS 


C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


PARAMETER ( NPRF=5 ) 

PARAMETER ( 2UNIT= ' FT ' , RUNIT=  *  NMI ' ) 

DIMENSION  N?T(11) ,RANGE(11) 

DIMENSION  2 (11, 30) ,R(11, 30) 

DIMENSION  21(11,310) ,RI(11, 310) 

ENTER  PROFILE  DATA 

INPUT  REQUIREMENTS :  1)  BOTH  RANGE  AND  HEIGHT  ARRAY  VALUES 

MUST  BE  IN  ASCENDING  ORDER. 

2)  FIRST  RANGE  ARRAY  ELEMENT  MUST  BE  ZERO. 

3)  ALL  HEIGHT  PROFILES  MUST  HAVE  THE  SAME 
VALUES  FOR  THE  FIRST  (AND  LAST)  ARRAY 
ELEMENTS . 

DATA  (N?T(I) ,1=1, 5) 728,28,28,24,24/ 

DATA  (RANGE (I) ,I=l,5)/0. , 80 . , 120 . , 160 . , 200 ./ 

DATA  (Z(1,J) ,J=l,28)/0. , 200 . , 300 . , 400 . , 500 . , 550 . ,600. ,650. ,700. , 

*  800. ,350. ,900. ,1000. , 12 00 . , 13 00 . , 14 00 , ,1500. ,1600. ,1800. , 

*  2000., 2200., 2400., 2600., 28 00., 3000., 3200., 34 00., 3500./ 

DATA  (R(1,J) , J=l, 28)/336. ,335., 334., 334., 332., 3 29., 326., 3 14., 

*  302. ,285. ,284. ,282. ,279. ,276. , 276 . , 276 . , 275 . , 273 . , 

*  272., 272., 271., 2 69., 2 67., 2 66., 266., 266., 2 64., 263./ 

DATA  (2 (2, J) , J=l,28)/0. , 200 . , 400 . , 500 . , 700. ,800. ,900. ,950. , 1000. , 

*  1050. , 1100. , 1200. , 1300. , 1400. , 1500. , 1550. , 1600. , 1700. , 1300. , 

*  2000. ,2200. ,2400. ,2600. ,2800. ,3000. ,3200. ,3400. ,3500./ 

DATA  (R(2,J) ,J=l,28)/337. ,335., 333., 332., 330., 324., 3 09., 304., 

*  296.  ,293. ,291. ,290. ,287. ,286. ,282. ,275. ,273. ,273. , 

*  272. ,272. ,271. ,269. ,267. ,266. ,266. ,266. ,264. ,263./ 

DATA  (2(3 ,J) , J=l,28)/0. , 200. ,300. ,400. , 600. , 700. ,800. , 1000. , 

*  1100. ,1200. ,1250. ,1300. ,1400. ,1450. ,1500. ,1550. ,1600. ,1800. , 

*  1900. ,2000. ,2100. ,2200. ,2400. ,2600. ,2800. ,3000. ,3300. ,3500./ 
DATA  (R(3,J) , J=l,28)/334. ,334. ,333. ,333. ,331. ,330. ,328. ,327. , 

*  3  27.  ,3  26. ,319. ,312. ,295. ,291. ,286. ,284. ,283. ,281. , 

*  276. ,271. ,270. ,269. ,269. ,269. , 266. , 265 . , 263 . , 263 ./ 

DATA  (2(4,J) ,J=l,24)/0. ,50. ,300. ,400. ,600. ,800. ,1000. ,1150. , 
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*  1200. ,14 00. ,1500. ,1600. ,1650. ,17 00. ,13 00. ,2000., 

*  22 00., 2400., 2500., 2 600., 28 00., 3000., 3300., 3500./ 
DATA  (R(4,J) , J=l, 24)/333 . ,332., 330., 3 29., 327., 325., 324 

*  320. ,305. ,297. ,292. ,289. ,287. ,282. ,273. , 

*  2 67., 271., 273., 272., 2 67., 2 65., 265., 261./ 

DATA  (Z (5, J) , J=l, 24)/0. ,2 00. ,3 00. ,400. ,600. ,800. ,1000. 

*  1400. , 1600. , 1700. , 1750. ,1850. ,2000. , 22 00 . , 2 4 00 . , 

*  2500. ,2600. ,2700. ,2800. ,2900. ,3100. ,3300. ,3500./ 
DATA  (R(5,J) ,J=1,24)/331. ,3 29., 327., 327., 324., 321., 320 

*  3 15., 3 07., 299., 301., 290., 280., 274., 270., 

*  273. ,274. ,273. , 2 65 . , 2 61 . , 2 64 . , 2 63 . , 2 62 . / 

OPEN  OUTPUT  FILE 

OPEN ( 9 , FILE= ' OUTFIT ' , STATUS= ' NEW ' ) 

CONVERT  N-UNITS  TO  M-UNITS 
DO  10  I  =  1,NPRF 
JM  =  NPT(I) 

DO  20  J  =  1,JM 

R(I,J)  =  R(I,J)  +  (2(1, J)*. 3048*. 157) 

20  CONTINUE 
10  CONTINUE 

CALL  RSMA  SUBROUTINE 

CALL  RSMA ( NPRF , ZUNIT , RUNIT , NPT , RANGE , Z , R , Z I , RI , I PMX ) 
WRITE  OUTPUT  FILE 
WRITE(9,30)  IPMX 

3  0  FORMAT  (' NUMBER  OF  LE"/ELS=  '  ,  14 ) 

DO  40  I  =  1,I?MX 

WRITE (9, ’ (10F7. 1) ' )  (ZI(M,I) ,RI(M, I) ,M=1,NPRF) 

40  CONTINUE 

STOP 


END  DEMO  PROGRAM 


END 


,  323 

1200 

,313 
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CCCCCCCCCCCCCCC  START  SUBROUTINE  RSMA  CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  RSMA ( NPRF , ZUNIT , RUNIT , NPT , RANGE , Z , R , Z I , RI , IPMX) 

CHARACTER*2  ZUNIT 
CHARACTER*3  RUNIT 

DIMENSION  NPT ( 11) , RANGE (11) 

DIMENSION  Z (11, 30)  ,R(11, 30) 

DIMENSION  ZI (11, 310) ,RI (11, 310) 

DIMENSION  DM(11,30) ,DDM(11,30) 

DIMENSION  IC(11, 30) ,MATCH(11,30) 

DIMENSION  ZM1(30) , ZM2 (30) 

DIMENSION  ZVAL(310) ,ZIF(11,310) ,ZIB(11,310) 

DIMENSION  ICAT1(36) ,ICAT2t36) 

DATA  ICATl/ 15, 14, 13, 12, 11,10,9,8,7,6, 5, 4, 3, 2, 1, 15, 15, 14, 

*  14,14,13,13,15,15,15,14,13,12,11,12,10,12,11,10,9,8/ 
DATA  ICAT2/ 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1,14, 13, 13, 

*  12,11,10,11,12,10,11,10,12,11,10,9,8,7,6,5,4,1/ 

DO  10  I  =  1,NPRF 
JM  =  NPT(I) 

CALCULATE  DM/DZ 
DM ( I , 1 )  =  0.0 
DO  30  J  =  2,JM 

DM(I,J)  =  (  R(I,J)  -  R(I,J-1)  )  /  (  Z(I,J)  -  Z(I,J-1)  ) 

IF(ZL^IT.EQ. ' FT' )  DM(I,J)  =  (  DM(I,J)  /  .3043  )  *  1000. 
IF(2UNIT.EQ. '  M')  DM(I,J)  =  DM(I,J)  *  1000. 

30  CONTINUE 

CALCULATE  CHANGE  IN  DM/DZ  AT  DATA  POINTS 

DDM(I,1)  =  0.0 
DDM(I,JM)  =0.0 

DO  40  J  =  2,JM-1 

DDM(I,J)  =  (  DM(I,J+1)  -  DM(I,J)  ) 


ono  o  o  n  non  no  non 


40  CONTINUE 


CATEGORIZE  DATA  POINTS  USING  DM/DZ  AND  D(DM/DZ) 

IC(I,1)  =  15 
IC(I,JM)  =  15 


DO  90  J  =  2,JM-1 


IF(DDM(I, J) .GT.-39. 
IF(DDM(I,J) .GT.-78. 
IF(DDM(I, J) .GE.78. ) 
IF(DDM(I,J) .GE.39. 


•AND-  DDM(I, J) .LT.39. )  GO  TO  100 
.AND.  DDM(I, J) .LE.-39. )  GO  TO  110 
GO  TO  120 

.AND.  DDM(I, J) .LT.78. )  GO  TO  130 


DDM(I,J)  <=  -78. 

IF(DM(I,J) .LT.157.)  GO  TO  91 
IF(DM(I,J-1) .LT.O. )  IC(I,J)  =-15 

IF(DM(I,J-1)  .GE.O.  .A-ND.  DM  ( I ,  J-^1 )  .  LT .  7  8  .  )  IC(I,J)  =-13 
IF(DM{I,J-rl)  .GE.78.  .AND.  DM  ( I ,  J+1)  .  LT .  157  .  )  IC(I,J)  =  -10 
IF(DM(I, J^l) .GE. 157. )  IC(I,J)  =-8 

91  IF(DM(I,J) .LT.78.  .OR.  DM ( I , J) . GE . 157 . )  GO  TO  92 
IF(DM(I,J+1) .LT.O. )  IC(r,J)  =-14 

IF(DM(I, J+1) .GE.O.  .AND.  DM ( I , J+1) . LT. 78 . )  IC(I,J)  =-11 

92  IF(DM(I,J) .LT.O.  .OR.  DM(I, J) .GE. 78 . )  GO  TO  93 
IF(DM(I,J+1) .LT.O.)  IC(I,J)  =-12 

93  IF(DM(I, J) .LT. 0. )  IC{I,J)  =  -9 
GO  TO  90 


DDM(I,J)  >  -73.  .AND.  DDM(I,J)  <=  -39. 

110  IF(DM(I, J) . LT. 157 . )  GOTO  111 
IF(DM(I, J^l) .GE. 73.  .AND.  DM(I, J+1) . LT. 157 . )  IC(I,J)  =-5 
IF(DM(I, J+1) .GE. 157 . )  IC(I,J)  =-l 

C 

111  IF(DM(I,J) .LT.78.  .OR.  DM(I, J) .GE. 157 . )  GO  TO  112 
IF(DM(I, J+1) .GE.O.  .AND.  DM ( I , J+1) . LT . 78 . )  IC(I,J)  =  -6 
IF(DM(I, J+1) .GE.78 .  .AND.  DM ( I , J+1) . LT . 157 . )  IC(I,J)  =  -2 

C 

112  IF(DM(I, J) .LT.O.  .OR.  DM(I, J) .GE. 78 . )  GO  TO  113 
IF(DM(I,J^1)  .LT.O.)  IC(I,J)  =  -7 

IF(DM(I, J+1) .GE.O.  .AND.  DM ( I , J+1) . LT . 78 . )  IC(I,J)  =-3 
C 

113  IF(DM(I, J) .LT.O. )  IC(I,J)  =-4 
GO  TO  90 

C 
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DDM(I,J)  >=  78. 

120  IF(DM(I,J+1) .LT.157.)  GO  TO  121 
IF(DM(I,J) .LT.O.)  IC(I,J)  =  15 

IF(DM(I, J) .GE. 0.  .AND.  DM ( I , J) . LT . 78 . )  IC(I,J)  =  13 
IF(DM(I, J) .GE.78.  .AND.  DM ( I , J) . LT. 157 . )  IC(I,J)  =  10 
IF(DM(I, J) .GE. 157. )  IC(I,J)  =  8 

121  IF(DM(I,J+1) .LT.78.  .OR.  DM ( I , J+1 ) . GE . 157 . )  GOTO  122 
IF(DM(I,J) .LT.O. )  IC(I,J)  =  14 

IF(DM(I, J) .GE.O.  .AND.  DM(I, J) . LT. 78 . )  IC(I,J)  =  11 

122  IF(DM(I,J-*-l)  .LT.O.  .OR.  DM  ( I ,  J+1)  .  GE .  78  . )  GOTO  123 
IF(DM(I, J) .LT.O. )  IC(I,J)  =  12 

123  IF(DM(I, J+1) .LT.O. )  IC(I,J)  =  9 
GO  TO  90 

DDM(I,J)  >=  39.  .AND.  DDM(I,J)  <  78. 

130  IF(DM(I, J+1) .LT.157.)  GO  TO  131 
IF(DM(I,J) .GE.78.  .AND.  DM ( I , J) . LT . 157 . )  IC(I,J)  =  5 
IF(DM(I, J) .GE. 157. )  IC(I,J)  *  1 

131  IF(DM(I,J+1) .LT.78.  .OR.  DM (I , J+1) . GE . 157 . )  GOTO  132 
IF(DM(I,J) .GE.O.  .AND.  DM ( I , J) . LT. 78 . )  IC(I,J)  =  6 
IF(DM(I,J) .GE.78.  .AND.  DM ( I , J) . LT. 157 . )  10(1, J)  =  2 

132  IF(DM(I, J+1) .LT.O.  .OR.  DM ( I , J+1) . GE . 78 . )  GOTO  133 
IF(DM(I,J) .LT.O.)  IC(I,J)  =7 

IF(DM(I, J) .GE. 0.  .AND.  DM ( I , J) . LT . 7 8 . )  IC(I,J)  =  3 

133  IF(DMCI, J+1) .LT.O. )  10(1, J)  =  4 
GO  TO  90 

100  10(1, J)  =  0 

90  oontinlt: 

10  OONTINUE 

MATCH  DATA  POINTS  OF  IDENTICAL,  THEN  SIMILAR,  CATEGORIES, 
BASED  ON  MATCH  HIERARCfff 

FTM  =  3.2808399 
C 

IF(ZUNIT.EQ.  'FT'  .AND.  RLTJIT.EQ.'  MI')  RFCT  =  5.28 
IF(ZUNIT.EQ. 'FT'  .AND.  RUNIT . EQ . ' NMI ' )  RFCT  =6.08 
IF(ZUNIT.EQ. 'FT'  .AND.  RUNIT. EQ.'  KM')  RFCT  =  3.23 


IZ 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


IF(ZUNIT.EQ. '  M'  .AND.  RUNIT.EQ.'  MI')  RFCT  =  1.61 
IF(ZUNIT.EQ. '  M'  .AND.  RUNIT. EQ . ' NMI ' )  RFCT  =  1.85 
IF(ZUNIT.EQ. '  M'  .AND.  RUNIT.EQ.'  KM')  RFCT  =1.00 

DO  200  M  -  1,NPRF-1 

RDEP  =  RFCT  *  (RANGE (M+1)  -  RANGE (M)) 

NMXl  =  NPT(M) 

NMX2  =  NPT(M+1) 

DO  210  ICAT  =  1,36 

DO  220  IP  =  1,4 

IF(IP.EQ.l)  J1  =  -1  *  ICATl(ICAT) 

IF(IP.EQ.l)  J2  =  -1  *  ICAT2(ICAT) 

IF(IP.EQ.2)  J1  =  ICATl(ICAT) 

IF(IP.EQ.2)  J2  =  ICAT2(ICAT) 

IF(IP.GE.3  .AND.  ICAT. LE. 15)  GO  TO  220 
IF(IP.EQ.3)  J1  =  -1  *  ICAT2(ICAT) 

IF(IP.EQ.3)  J2  =  -1  *  ICATl(ICAT) 

IF(I?.EQ.4)  J1  =  ICAT2(ICAT) 

IF(IP.EQ.4)  J2  =  ICAT1{ICAT) 

DO  230  N1  =  1,NMX1 

IF(IP.EQ.l  .AND.  ICAT.EQ.l)  MATCH(M,N1)  =  0 

IF(IC(M,N1) .NE.Jl)  GO  TO  230 
IF(MATCH(M,N1) .NE.O)  GO  TO  230 

ICNT  =  0 

IF(ZUNIT.EQ. '  M')  GO  TO  233 

IF(Z (M,N1) .LT. 100. )  WNDW  =  RDEP  +  200. 

IF(Z(M,N1) .GE.IOO.)  WNDW  =  RDEP  + 

*  (100.*(ALOG10(Z(M,N1) ) ) ) 

GO  TO  237 

233  IF(Z(M,N1) .LT. (lOO./FTM) )  WNDW  =  RDEP  +  (200. /FTM) 

IF(Z (M,N1) .GE. (lOO./FTM) )  WNDW  =  RDEP  + 

*  ( (lOO./FTM) *(ALOG10( (2(M,N1) *FTM) ) ) ) 

237  ZMIN  =  Z(M,N1)  -  WNDW 

ZMIN  =  MAX(Z(M,1) ,ZMIN) 

ZMAX  =  Z(M,N1)  +  WNDW 
ZMAX  =  MIN (Z(M, NMXl) ; ZMAX) 
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oooo  n  o  n  ooo 


c 

c 

c 


c 


c 


DO  240  N2  =  1,NMX2 

IF(Z (M+1,N2) .LT.ZMIN  .OR.  Z (M+1 , N2 ) . GT . ZMAX)  GO  TO  240 

IF(IC(M+1,N2) .NE.J2)  GO  TO  240 

ICNT  =  ICNT  +  1 
IF(ICNT.GT. 1)  GO  TO  250 

MATCH (M,N1)  =  N2 

CALL  MTHCHK ( MATCH, M,NMX 1,1 FLAG) 

IF  (IFLAG.EQ.O)  GO  TO  260 

IFLAG  =  1  REJECT  MATCH;  RESET  COUNTER  AND  MATCH  VALUE 


ICNT  =  0 
MATCH (M,N1)  =  0 
GO  TO  240 

260  MTHINT  =  N2 

ZDIF  =  ABS(2(M,N1)  -  2(M+1,N2)) 

GO  TO  240 

250  ZDIFN  =  ABS(Z(M,N1)  -  2(M+1,N2)) 

IF(2DIFN.GE.ZDIF)  GO  TO  240 

MATCH{M,N1)  =  N2 
CALL  MTHCHK (MATCH, M,NMX1, IFLAG) 

IF  (IFLAG.EQ.O)  GO  TO  270 

IFLAG  =  1  REJECT  MATCH;  RESET  COUNTER  AND  MATCH  VALUE 

ICNT  =  ICNT  -  1 
MATCH (M,N1)  =  MTHINT 
GO  TO  240 
C 

270  MTHINT  *  N2 

ZDIF  =  ZDIFN 
C 

240  CONTINUE 

C 

230  CONTINUE 

C 

220  CONTINUE 

210  CONTINUE 

C 

200  CONTINUE 
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o  n  o  o  n 


PERFORM  FORWARD,  THEN  BACKWARD,  MATCHING;  ON  COMPLETION, 
ALL  PROFILE  DATA  ARE  MATCHED  SUCH  THAT  ALL  PROFILES  HAVE 
SAME  NUMBER  OF  POINTS. 


C 


C 

c 


c 


c 


DO  400  ID  =  1,2 


IF(ID.EQ.l)  Ml 
IF(ID.EQ.l)  M2 
IF(ID.EQ.l)  M3 
IF(ID.EQ.2)  Ml 
IF(ID.EQ.2)  M2 
IF(ID.EQ.2)  M3 


1 

NPRF-1 

1 

NPRF 

2 


-1 


DO  405  M  =  Ml, M2, M3 


IF(ID.EQ.l)  NMX  =  NPT(M) 
IF(ID.EQ.2)  NMX  =  NPT(M-l) 
NMP  =  0 


DO  410  N  =  1,NMX 

IF(ID.EQ.2)  GO  TO  415 

L  =  MATCH (M,N) 

IF(MATCH(M,N) .EQ.O)  GO  TO  410 
NMP=  NMP  +  1 
2M1(NMP)  =  Z(M,N) 

2M2(NMP)  =  Z(M+1,L) 

GO  TO  410 


C 

415 


C 

410 

C 


C 


c 


L  =  MATCH (M-1,N) 

IF(MATCH(M-1,N) .EQ.O)  GO  TO  410 
NMP  =  NM?  +  1 
ZMl(NMP)  =  Z(M,L) 

ZM2(NMP)  =  Z(M-1,N) 

CONTINUE 

IF(M.EQ.l  .OR.  M.EQ.  NPRF)  NPI  =  NPT(M) 
IF(M.GT.l  .AND.  M.LT.NPRF)  NPI  =  IP 
NMP  =  0 
NIP  =  0 

DO  420  N  =  1,NPI-1 

IF(M.EQ.l  .OR.  M.EQ. NPRF)  ZVAL(N)  =  Z(M,N) 
IF(N.EQ.l)  GO  TO  430 
IF(ZVAL(N) .EQ.ZT)  GO  TO  430 

NIP  =  NIP  +  1 
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* 


it 


c 

430 


C 

420 

C 


c 


c 


IF(ID.EQ. 1) 

ZIF(M+1,NIP)  =  ( ( (ZVAL(N)-ZB)/ZDIF1) *ZDIF2) 
IF(ID.EQ. 2) 

ZIB(M-1,NIP)  =  ( ( {ZVAL(N)-ZB)/ZDIF1) *ZDIF2) 
GO  TO  420 


NMP  =  NMP  +  1 
ZB  =  ZMl(NMP) 

ZT  =  ZM1(NMP+1) 

ZDIFl  =  ZM1(NMP+1)  -  ZM1{NMP) 
ZDIF2  =  ZM2(NMP+1)  -  ZM2 (NMP) 


CONTINUE 
IP  =  1 

IF(ID.EQ.l)  ZVAL(l)  =Z(M+1,1) 
IF(ID.EQ.2)  ZVAL(l)  =  Z(M-1,1) 


IF(ID.EQ.l)  NMX  =  NPT(M+1) 
IF(ID.EQ.2)  NMX  =  NPT(M-l) 
NINT  =  1 


DO  440  N  =  2, NMX 


C 

C 


c 


c 

450 

C 

C 

460 


C 

440 

C 


IFLAG  =»  0 
NUMINT  =  0 

IF (NINT. GT. NIP)  GO  TO  460 
DO  450  L  =  NINT, NIP 

IF(ID.EQ.l  .AND.  ZIF (M+1 , L) . GE . Z (M+1 , N) ) 
IF(ID.EQ.2  .AND.  ZIB (M-1 , L) . GE . Z (M-1 , N) ) 

IFLAG  =  1 

NUMINT  =  NUMINT  +  1 
IP  =  IP  -t-  1 

IF(ID.EQ.l)  ZVAL(IP)  =  ZIF(M+1,L) 
IF(ID.EQ.2)  ZVAL(IP)  =  ZIB(M-1,L) 

CONTINUE 

IF  ( IFLAG.  EQ.  1)  NINT  =  NINT  -t-  NUMINT 
IP  =  IP  +  1 

IF(ID.EQ.l)  ZVAL(IP)  =  Z(M+1,N) 

IF(ID.EQ.2)  ZVAL(IP)  =  Z(M-1,N) 

CONTINUE 


+  ZM2(NMP) 
+  ZM2(NMP) 


GO  TO  450 
GO  TO  450 
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c 


c 


c 

c 


c 

c 


c 

c 


c 

c 


c 

c 

c 


c 

c 


c 


IF(ID.EQ.2)  GO  TO  470 

IF(M.NE. (NPRF-1) )  GO  TO  405 
IPMX  =  IP 

DO  480  N  =  1,IP 

ZI(M+1,N)  =  ZVAL(N) 

480  CONTINUE 

GO  TO  405 

470  IF(M.NE.2)  GO  TO  490 

DO  500  N  =  1,IP 

ZI(M-1,N)  =  ZVAL(N) 

500  CONTINUE 

GO  TO  405 

490  NZIF  =  IPMX  -  IP 

1  =  1 

ZI(M-1,1)  =  Z(M-1,1) 

NINT  =  1 

DO  510  N  =  2, IP 

IFLAG  =  0 
NUMINT  =  0 

IF(NINT.GT.NZIF)  GO  TO  520 

DO  530  L  =  NINT, NZIF 

IF(ZIF(M-1,L) .GE.ZVAL(N) )  GO  TO  530 
IFLAG  =  1 

NUMINT  =  NUMINT  +  1 
1  =  1  +  1 

ZI(M-1,I)  =  ZIF(M-1,L) 

530  CONTINUE 

IF ( IFLAG. EQ. 1)  NINT  =  NINT  +  NUUINT 

520  1=1+1 

ZI(M-1,I)  =  ZVAL(N) 

510  CONTINUE 

405  CONTINUE 
400  CONTINUE 
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ooo  o  oooooo  noo  no 


REMOVE  REPETITIONS  IN  ZI  ARRAYS . 

DO  540  M  =  1,NPRF 

IMX  =  IPMX  -  1 

DO  550  I  =  1,IMX 

IF(I.GT.IMX)  GO  TO  550 
IF(2I(M,I+1) .NE.ZI(M,I) )  GO  TO  550 

II  =  I  +  1 

DO  560  J  =  II, IMX 

ZI(M,J)  =  ZI(M,J+1) 

560  CONTINUE 

IMX  =  IMX  -  1 

550  CONTINUE 
540  CONTINUE 

IPMX  =  IMX  +  1 

CALCULATE  RI  VALUES  FOR  ZI  POINTS,  USING  LINEAR  INTERPOLATION. 

DO  600  M  =  1,NPRF 
NMX  =  NPT(M) 

INIT  =  1 

DO  610  N  =  1,NMX-1 
NUm.NT  =  0 

RB  =  R(M,N) 

RT  =  R(M,N+1) 

RDIF  «  RT  -  RB 
ZB  =  Z(M,N) 

ZT  »  2(M,N+1) 

ZDIF  =  ZT  -  ZB 

DO  620  I  *  INIT, IPMX 

IF(2I(M,I) .EQ.ZT)  GO  TO  630 

RI(M,I)  =  RB  (  (  (ZI(M,I)-ZB)/ZDIF)  *RDIF) 

NUMINT  »  NUMINT  -f  1 
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620  CONTINUE 

C 

C30  INIT  =  INIT  +  NUMINT 

IF(I.EQ. IPMX)  =  RT 

C 

610  CONTINUE 

600  CONTINUE 
C 

RETURN 

END 

C 

SUBROUTINE  MTHCHX (MATCH , M, N , IFLAG) 
C 

DIMENSION  MATCH(11,30) 

C 

IFLAG  =  0 
C 

LOW  =  MATCH(M, 1) 

C 

DO  100  I  =  2,N 
C 

IF(MATCH(M, I) .EQ.O)  GO  TO  100 
C 

NEXT  =  MATCH (M, I) 

IF (NEXT. LE. LOW)  GO  TO  200 
LOW  *  NEXT 
GO  TO  100 
C 

200  IFLAG  =  1 

C 

100  CONTINUE 

C 

RETURN 

END 
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OUTPUT  -  DEMO  PROGRAM 


NUMBER  OF  LEVELS*  98 


.0 

336.0 

.0 

337.0 

.0 

334.0 

.0 

333 . 0 

.0 

331.0 

14.5 

336.6 

20.3 

337.8 

34.8 

335.7 

33 . 3 

333.9 

200.0 

338.6 

21.7 

336.9 

30.4 

338.2 

52.2 

336.5 

50.0 

334.4 

300.0 

341.4 

58.5 

338.5 

81.9 

340.1 

140.5 

340.7 

134.6 

337.8 

400.0 

346.1 

83.3 

339.6 

116.7 

341.4 

200.0 

343.6 

191.7 

340.0 

467.4 

348.4 

125.0 

341.4 

175.0 

343.6 

300.0 

347.4 

287.5 

343.9 

580.7 

352.1 

130.4 

341.6 

182.6 

343.9 

313.0 

348.0 

300.0 

344.4 

595.5 

352.6 

132.1 

341.7 

184.9 

344.0 

317.1 

348.2 

303.8 

344.5 

600.0 

352.7 

142.9 

342.1 

200.0 

344.6 

342.9 

349.4 

328.6 

345.4 

629.2 

353.7 

166.7 

343 . 1 

233 . 3 

345.8 

400.0 

352.1 

383.3 

347.5 

693.9 

355.8 

173.9 

343.5 

243.5 

346.2 

417.4 

352.8 

400.0 

348.1 

713 . 6 

356.4 

200.0 

344.6 

280.0 

347.6 

480.0 

355.2 

460.0 

350.4 

784.5 

358.3 

205.7 

344.8 

288.0 

347.9 

493.6 

355.7 

473 . 1 

350.9 

800.0 

359.3 

250.0 

346.5 

350.0 

350.2 

600.0 

359.7 

575.0 

354.3 

920.5 

364.4 

260.9 

346.9 

365.2 

350.8 

626.1 

360.7 

600.0 

355.7 

950.0 

365.7 

279.3 

347.6 

391.0 

351.8 

670.2 

362.4 

642.3 

357.3 

1000.0 

367.9 

285.7 

347.8 

400.0 

352.1 

685.7 

363.0 

657.1 

357.9 

1017.5 

363.5 

291.7 

348.0 

408.3 

352.5 

700.0 

363.5 

670.8 

353.4 

1033.7 

369.1 

300.0 

348.4 

420.0 

352.9 

720.0 

364.1 

690.0 

359.1 

1056.4 

370.0 

333.3 

350.0 

466.7 

354.7 

800.0 

366.3 

766.7 

362.0 

1147.0 

373.4 

347.8 

350.6 

487.0 

355.4 

834.8 

367.8 

800.0 

363.3 

1186.4 

374.9 

352.8 

350.9 

494.0 

355.7 

846.8 

368.3 

811.5 

363.8 

1200.0 

375.4 

357.1 

351.1 

50C.0 

355.9 

857.1 

368.7 

821.4 

364.2 

1211.7 

375.3 

400.0 

353 . 1 

560.0 

358.2 

960.0 

373.1 

920.0 

368.4 

1323.2 

379.6 

416.7 

353 . 6 

583.3 

359.1 

1000.0 

374.9 

958.3 

370.1 

1373 . 5 

381.1 

426.4 

353.9 

597.0 

359.6 

1023.4 

376.0 

980.8 

371.0 

1400.0 

332.0 

434.8 

354.1 

608.7 

360.0 

1043.5 

376.9 

1000.0 

371.9 

1422.7 

382.2 

458.3 

354 . 3 

641.7 

361.3 

1100.0 

379.6 

1054.2 

374.1 

1486.7 

382.7 

500.0 

355.9 

700.0 

363.5 

1200.0 

383.4 

1150.0 

378.0 

1600.0 

383 . 6 

550.0 

355.3 

750.0 

362.9 

1250.0 

378.8 

1175.0 

377.7 

1617.9 

383.0 

600.0 

354.7 

800.0 

362.3 

1300.0 

374.2 

1200.0 

377.4 

1635.7 

382.4 

630.0 

348.9 

830.0 

359.2 

1330.0 

370.5 

1290.0 

375.0 

1700.0 

380.4 

650.0 

345.1 

850.0 

357.2 

1350.0 

368.1 

1350.0 

373.4 

1742.9 

384.1 

653.3 

344.5 

853.3 

356.8 

1353.3 

367.7 

1360.0 

373.1 

1750.0 

384.7 

666.7 

341.9 

866.7 

355.5 

1366.7 

366.1 

1400.0 

372.0 

1778.6 

383.0 

700.0 

335.5 

900.0 

352.1 

1400.0 

362.0 

1500.0 

368.8 

1850.0 

378.5 

725.0 

332.4 

950.0 

349.5 

1450.0 

360.4 

1575.0 

368.6 

1925.0 

377.1 

750.0 

329.4 

1000.0 

343.9 

1454.2 

360.2 

1581.3 

368.6 

1931.3 

377.0 

775.0 

326.3 

1050.0 

343.2 

1458.3 

360.0 

1587.5 

368.6 

1937.5 

376.9 

800.0 

323.3 

1100.0 

343.6 

1462.5 

359.7 

1593.8 

368.6 

1943.8 

376.8 

825.0 

324.0 

1150.0 

345.5 

1466.7 

359.5 

1600.0 

368.6 

1950.0 

376.6 

850.0 

324.7 

1200.0 

347.4 

1470.8 

359.3 

1606.3 

368.5 

1956.3 

376.5 

900.0 

325.1 

1233.3 

348.0 

1473.6 

359.2 

1610.4 

368.4 

1960.4 

376.5 
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1000 . 0 

326.9 

1300.0 

349 . 2 

1479.2 

353.9 

1618.8 

368.3 

1968 . 3 

376.3 

1113 . 6 

330.6 

1400.0 

353 . 0 

1487 . 5 

358.4 

1631.3 

36S  .  2 

1931 . 3 

376.1 

1200.0 

333.4 

1476.0 

353.6 

1493.8 

358.1 

1640.8 

368 , 1 

1990.8 

375 . 9 

1227 . 3 

334.7 

1500.0 

353 . 8 

1495.8 

358.0 

1643 . 8 

363.0 

1993 . 3 

375.8 

1234 . 1 

337.4 

1550.0 

349 . 2 

1500.0 

357.8 

1650.0 

368.0 

2000 . 0 

375.7 

1289.3 

337.7 

1554 . 5 

349 . 2 

1545.5 

358.1 

1700 . 0 

363.4 

2036.4 

376.4 

1289.8 

337.7 

1555.0 

349.2 

1550.0 

353.2 

1705.0 

368 . 3 

2040 . 0 

376.4 

1295.5 

338 . 0 

1560.0 

349.3 

1600.0 

359.6 

1760.0 

368.2 

2080.0 

377 . 1 

1299.6 

338 . 2 

1563 . 6 

349.3 

1636.4 

360.9 

1300.0 

363.1 

2109 . 1 

377.7 

1300.0 

338 . 2 

1564 . 0 

349.3 

1640.0 

361.1 

1804 . 0 

363.1 

2112 . 0 

377.7 

1312 . 5 

338.8 

1575.0 

349.4 

1750 . 0 

365.2 

1925.0 

363.5 

2200.0 

379 . 3 

1318.2 

339.1 

1580.0 

349 . 4 

1800.0 

367.1 

1930 . 0 

368.7 

2240 . 0 

380 . 4 

1320.2 

339.2 

1581.8 

349 . 4 

1813.2 

367.1 

2000.0 

363.7 

2254 . 5 

380.8 

1329.5 

339 . 6 

1590.0 

349.5 

1900.0 

366.9 

2090.0 

370.3 

2320.0 

382 . 6 

1340.9 

340.2 

1600 . 0 

349 . 6 

2000.0 

366.7 

2200 . 0 

372 . 3 

2400.0 

384.3 

1400.0 

343 . 0 

1652 . 0 

352 . 1 

2041.1 

363.3 

2235.6 

374.7 

2435.6 

387 . 6 

1454 . 5 

345.1 

1700.0 

354.4 

2078.9 

369.7 

2268.4 

376.9 

2463.4 

390.2 

1484 . 8 

346.2 

1726.7 

3  55.4 

2100.0 

370.5 

2286.7 

378 . 2 

2486.7 

391.6 

1500.0 

346.8 

1740.0 

355.9 

2110.5 

370.9 

2295.3 

373 . 3 

2495 . 8 

392 . 3 

1507.0 

347 . 0 

1746.2 

356.1 

2115.4 

371.1 

2300.0 

379 . 1 

2500.0 

392 . 6 

1568.2 

343 . 7 

1300 . 0 

358.1 

2157.9 

372.7 

2336.8 

381.6 

2536.3 

394 . 3 

1600.0 

349.6 

1828.0 

359.5 

2130.0 

373 . 5 

2356.0 

332.9 

2556.0 

395.9 

1623.8 

350.3 

1853 . 3 

360.7 

2200,0 

374.3 

2373 . 3 

334 . 0 

2573 . 3 

396.9 

1673 . 1 

352.7 

1892 . 3 

362 . 6 

2230.8 

375.3 

2400.0 

385.8 

2600 . 0 

398.4 

1795.5 

357.9 

2000.0 

367.7 

2315.3 

379.8 

2473.7 

390 . 8 

2673.7 

404.9 

1800.0 

358.1 

2004 . 0 

367.9 

2318.9 

330.0 

2476.4 

391.0 

2676.4 

405.1 

1339.2 

360.0 

2038 . 5 

369 . 4 

2346.2 

381.3 

2500.0 

392.6 

2700.0 

407.2 

1916.7 

363.7 

2106.7 

372 . 3 

2400.0 

383.8 

2546.7 

394.4 

2737.3 

404 . 1 

2000.0 

367.7 

2180.0 

375.4 

2457.9 

386.6 

2596.8 

396.3 

2777.5 

400.8 

2005.2 

367.9 

2184.6 

375.6 

2461.5 

386.8 

2600.0 

396.4 

2780.0 

400.6 

2022.7 

363.7 

2200.0 

376.3 

2473.7 

337.4 

2610.5 

396.7 

2788 . 4 

399.9 

2046.8 

369.7 

2221.2 

377,1 

2490.4 

338.2 

2525.0 
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